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ABSTRACT 

The applicability of Richardson's "four-thirds law" to turbulent 
horizontal diffusion, the dependence of diffusion on the weight of the 
diffusers, and the effects of varying the rejection level of the data 
were investigated. Diffusers of identical size and shape but weighing 
13 or 15 pounds were photographed from a Navy US2-A aircraft. The 
data were collected in Monterey Bay on scales from 25 meters to 200 
meters, with water depths ranging from 220 to 270 feet. 

The results indicate a nearly constant value of k (the constant 

4/3 

in Richardson's "four-thirds law", ie . F(l) = kl ) for Monterey Bay, 
although in nearly all cases the slope of the best-fit line was 
greater than predicted by Richardson (1926). The weight effect 
remains unsettled. Varying the rejection criterion has definite 
effects on k and on the slope of the best-fit line. 
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I . INTRODUCTION 



A, HORIZONTAL TURBULENT DIFFUSION AND RICHARDSON’S "4/3" LAW 

Richardson (1926), after analysis of atmospheric experimental data, 
noted that in molecular diffusion the motion of each molecule is inde- 
pendent of its immediate neighbors, whereas, in turbulent diffusion it 



Sv 



a 2 v 



is not. Then the classical Fickian equation = K — where v is 

Sx 

the concentration of a diffusing substance, x is position, t is time, 
and K is the diffusivity (a measure of the rate of diffusion), cannot 
be applied to turbulent diffusion, for it predicts a constant rate of 
diffusion with x. Stommel (1949) pointed out that using the Fickian 
equation the probability of two particles initially a distance of 
apart being a distance of l apart at some later time, would be 

p(, °- li> “’(.--e-;- 

© 

Therefore, the Fickian equation leads to the result that the probability 

of a pair of particles a distance of 1^ apart being a distance of 

2 

apart after an interval of time t depends upon (l^-l^) only, and not 
upon either 1^ or l^. This -is clearl)T in direct contradiction to the 
atmospheric observations. There seems to be no way in which the Fickian 
equation can be modified to overcome this discrepancy. 

On the basis of atmospheric diffusion experiments, Richardson (1926) 
postulated the equation 



a 

a? = a? 



F <*> IS- 



Here, l is the distance between particles and is called the "neighbor 
separation". The number of particles which have neighbors with neighbor 
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separations of between l and l + dl is q ( l ) d t . The quantity q ( t ) may 
then be called "neighbor concentration". The term F(l) is the "neighbor 
diffusivity" and is analogous to the diffusivity (K) in the Fickian 
equation . 

Richardson empirically determined that neighbor diffusivity was 

4/3 

related to neighbor separation by F(l) = kl where k is a constant. 
This is commonly refereed to as the "four-thirds law". Stommel (1949), 
on the basis of data obtained from aerial photography, suggested that 
the four-thirds law was also applicable to oceanic diffusion. Then, 
assuming the turbulence causing the diffusion is within the inertial 
subrange, it can be shown that k = CE 9 where C is a universal 
dimensionless constant and E is the turbulent energy dissipation per 
unit mass . 

Many individual studies have been conducted to test Richardson's 

four-thirds law. While general agreement exists as to the validity of 

the law, the data have not been sufficient to actually confirm it. 

Attempts have been made to determine k , but very little work has been 

done on evaluating C. The reported values of k range from 0.003 to 
2/3 -1 

0.09 cm sec . (These values were found by passing a best-fit line 
of 4/3 slope through the data points plotted on logarithmic paper.) 

B. A METHOD FOR THE DETERMINATION OF THE NEIGHBOR DIFFUSIVITY AS A 
FUNCTION OF THE NEIGHBOR SEPARATION 

Consider two particles initially separated by a distance l at 

time t^. At a later time t. the separation of particles is L. If 

l -l 1 

the time interval t^t^ is chosen so that the ratio — is small, 
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Richardsons equation can be written as 



dt 



FO 0 ) 




which has the solution 



qOp 



cons t . 




4t E ( 1 0 ) J * 



Thus the neighbor separation has a Gaussian distribution. The standard 

2 

deviation of from the mean is o = JlLt F (t . Thus F(Iq) = *2At' 

- N 2 — 

The best estimate of a results in F(t)= £ (At ) where l is an average 

(N- 1) 2At 

of l and (since all pairs of particles would not have the same 
initial separations), and N is the number of pairs averaged over. 



C. OBJECTIVES OF THE STUDY 

This study was conducted with three main objectives in mind. The 
first was to test the applicability of Richardson’s four-thirds law to 
turbulent horizontal diffusion in the ocean under various field con- 
ditions. An investigation of the dependence of diffusion on the weight 
of the particles was the second objective. This was studied by com- 
paring the diffusion of 13 pound and 13 pound plywood boards. (Earlier 
work had lead to the belief that even such a slight difference might be 
important.) The third objective was to study the effects of varying the 
rejection level of the data. That is, in most of the previous work the 

data were considered acceptable if the relative change in separation 
2 Li 

R = " — — was on the order of 0.10. By giving R various upper limits 

0 1 

it was felt that the data would be more useful and the effects, if any, 
of the different limits could be determined. This was accomplished by 
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calculating the neighbor diffusivity and the neighbor separation for 
all pairs with a relative change in separation of less than 0.05, less 
than 0.10, and less than 0.15. 
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II. EQUIPMENT AND PROCEDURES 



A. TECHNIQUES USED IN COLLECTION OF DATA 

The data used in this study were collected by Champion (1968)* The 
particles were 4* x 4 f x l/8 n plywood sheets, which were treated with 
sealer and marked for easy identification. The weights of the sheets 
were altered with the use of lead strips to within 0.1 pound of the 
prescribed weights of 13 and 15 pounds. The 13 pound sheets (A) were 
marked with a large red dot on a white background. The 15 pound 
sheets (B) were marked with two red horizontal stripes on a white 
background . 

The sheets were photographed with a T-ll Fairchild aircraft mapping 
camera installed in the rear section of a Navy US2-A aircraft (Philipps, 
1968). The camera was equipped with a 6-inch f:6.3 (Bausch and Lomb) 
class T metrogon lens that produced 9-1/2" x 10-1/4" negatives. The 
altitude of the camera to the nearest foot and time to the nearest 
second were automatically recorded on each negative. 

The data were collected about 2300 yards offshore in Monterey Bay, 
with water depths ranging from 220 feet to 270 feet. The approximate 
initial location of the patterns was 36°-38.4 , N and 121°-53.0 , W. 

Circular patterns of about 280 feet in diameter were used. The time 
between each pass of the aircraft over the pattern was approximately 
two minutes. During each pass a series of three to six photos were 
taken to improve the possibility of getting the entire pattern into 
at least one photo for that pass. If a specific sheet was missing 
from a certain pass, it was located in the next pass and the analysis 
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continued. The final shape of the patterns for both runs was approxi- 
mately elliptical with the major axis perpendicular to the direction of 
the mean current. 
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III. DATA ANALYSIS 



Appendix A contains the pertinent data for the two flights used in 
this study. Data were analyzed from 24 September 1968 (Run I) and 
1 October 1968 (Run II). The temperature structure of the upper layer 
on these days is shown in Figure 1. The negatives were used as 
transparencies and projected with a Travel-Graph overhead projector 
on an x-y grid. The grid was a sheet of 5 f x 5 1 plywood painted white, 
with coordinates marked by telephone cable wire spaced at 6-cm. inter- 
vals. After the best transparencies were selected for each pass, each 
sheet was properly located and marked. The transparencies were then 
projected on the x-y grid. (The enlargement factor was 5.5.) The 
x-y coordinates of each sheet were determined, and a digital computer 
was used to calculate the initial separations Oq)* After removal of 
all duplicate data cards, another program computed the mean separations 



l ^+1 . 



(^— 2 — ^)j the change in separation (At), the relative change in 
separation (R) , and neighbor diffusivity J for all pairs. 

Distances on the negative were converted to actual distances using 
the relationship: 



IMAGE = GROUND COVERED . 
FOCAL LENGTH ALTITUDE 



Histograms of the mean separations were drawn using three rejection 
factors: (i) R < 0.05, (ii) R < 0.10, and (iii) R < 0.15. Class 

intervals of the mean separations were established and the average 
value (over the class interval) of the mean separation and neighbor 



r N 2i 

diffusivity £ (At ) were found using a desk calculator. Very few 
(N-l)2/\t 
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values were added by changing the rejection factor from 0.10 to 0.15. 
Thus, only data with (i) R < 0.05 and (ii) R < 0.10 were considered. 



18 



TEMPERATURE (°C) 





W 






▼H 








w 








Dd 

Dd 


jg 


• 

w 


£ 


O 




> 


w 


M 


p-. 




Q 




s 

&H 







<ao 




19 



25 SEPT 
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IV. RESULTS AND CONCLUSIONS 



A. CONSISTENT RESULTS OF THE DATA 

A total of 6,399 observations were used, covering scales from 25 
meters to 200 meters. Appendix B gives the results obtained from these 
data. Eighty percent confidence limits were calculated for F(l) using 
the number of values averaged over minus one as the degrees of freedom 
of a chi-square distribution. The calculated values of F(l) and l were 
plotted on logarithmic paper. A best-fit line, in the least squares 
sense, and a best-fit line of four-thirds slope were drawn through the 
data points (Figures 2-9). 

Comparing the results of sheets A and B from Run I with those from 
Run II clearly demonstrates the consistent results of the data col- 
lected on the two different days. The best-fit lines for sheets A 
nearly coincide (Figure 10), and the best-fit lines for sheets B, 
even though they are separated by a small vertical distance, have 
approximately the same slope (Figure 11). The wave conditions in the 
two runs were quite different (see below). 

In nearly all cases the slope of the best-fit line was greater than 
predicted by Richardson (1926). Table I compares these results with 
Snyder (1967) and Philipps (1968), for a common rejection factor of 
0.10. Three of the four slope values determined from Runs I and II 
were larger than the value obtained by Philipps. A somewhat larger 
difference in slopes is noted when compared with Snyder. These 
differences will be discussed in a later section. 
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TABLE I 



SLOPE VALUES FOR THE BEST-FIT 
(FOR A REJECTION FACTOR OF 0. 


LINE 

,10) 


SOURCE 


SLOPE 


Snyder (March, 1967) Paper sheets 


1.46 


Philipps (April, 1968) 10 lb. plywood sheets 


1.612 


Run I 13 lbs . 


1.62 


Run I 15 lbs . 


1.27 


Run II 13 lbs. 


1.86 


Run II 15 lbs. 


1.71 


The values of k determined by passing a best- 

\ 


-fit line of four- 


thirds slope through the data are shown in Table 


II (using a common 



rejection factor of 0.10). The near equality of the k values for 



Monterey Bay determined by various investigators 
little with time. 


suggest that k varies 


TABLE II 





k VALUES FOR THE BEST-FIT LINE OF FOUR-THIRDS SLOPE 
(FOR A REJECTION FACTOR OF 0.10) 



SOURCE 


v ^ 2/3 -lx 

k (cm sec ) 


Snyder (March, 1967) Paper sheets 


0.006 



Philipps (April, 1968) 10 lb. plywood sheets 0.006 + 0.002 



Run I 13 lbs . 


0.004 


Run I 15 lbs. 


0.003 


Run II 13 lbs . 


0.005 


Run II 15 lbs. 


0.005 
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FIGURE 2 

NEIGHBOR DIFFUSIVITY F 
vs. — 

NEIGHBOR SEPARATION J- 

(Run It 13 sheets » R^0«05) 
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FIGURE 3 

NEIGHBOR DIFFUSIVITY F 
vs. _ 

NEIGHBOR SEPARATION JL 

(Run I, 13 lb. sheets, R<0.10) 



o 




2x1 0^ X (cm) 10 4 
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FIGURE 4 

NEIGHBOR DIFFUSIVITY F 
vs. _ 

NEIGHBOR SEPARATION JL 

(Run I, 15 lb. sheets, R<0.05) 
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FIGURE 5 

NEIGHBOR DIFFUSIVITY F 
vs. _ 

NEIGHBOR SEPARATION JL 

(Run I, 15 lb. sheets, R^D.IO) 
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FIGURE 8 

NEIGHBOR DIFFUSIVITY F 
vs. _ 

NEIGHBOR SEPARATION JL 

(Run II, 15 lb. sheets, R<0.05) 
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FIGURE 9 

NEIGHBOR DIFFUSIVITY F 
vs. _ 

NEIGHBOR SEPARATION J- 

(Run II, 1 5 lb. sheets, R<0.10) 
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FIGURE 10 

NEIGHBOR DIFFUSIVITY F 
vs. _ 

NEIGHBOR SEPARATION X 

(Run I, 13 lb. sheets, R<0.05) A 

(Run II, 13 lb. sheets, R<0.05) O 
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FIGURE 11 

NEIGHBOR DIFFUSIVITY F 
vs. j- 

NEIGHBOR SEPARATION X 

(Run I, 15 lb. sheets, R<0.05) O 

(Run II, 15 lb. sheets, R<t).05) A 
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B. WEIGHT EFFECT 



As has been pointed out by other authors, eddies with length 
scales less than the size of the dispersing particles cannot effect 
the dispersion (eg., Ozmidov, 1957). There is another side to this. 
Due to inertial effects, heavier particles will be less sensitive to 
rapid turbulent fluctuations than light particles. If the turbulent 
time and space scales are roughly proportional (ie., smaller eddies 
wiggle faster), then inertial effects will be more pronounced at 
smaller scales, and F(l) at these scales for heavy particles should 
be less than F ( t ) for light particles. As the scales increase, the 
effect should diminish. Thus the slope of the F, l curve should 
increase with weight. 

The effect of size was investigated by Ozmidov (1957), who con- 
firmed the speculation given above (although he did not mention the 
weight effect). A comparison between the data of Philipps, Snyder, 
and that analyzed here shows similar effects (Figures 12 and 13), 

(The slope change for Philipps is probably due mostly to weight; 
that for Snyder is unclear.) In an effort to isolate the weight 
effect, 13 lb. and 15 lb. plywood sheets of identical size and shape 
were used, as described above. However, the effect (although slight) 
was opposite to that conjectured above (Table III). This suggests 
that time and space scales are not proportional, and we must regard 
the problem of weight effect as unsettled. 

C. EFFECTS OF THE REJECTION CRITERION 

There has been little consideration given to the effects of 
varying the rejection criterion on k and on the slope of the best-fit 
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FIGURE 12 

NEIGHBOR DIFFUSIVITY F 

vs. — O 

NEIGHBOR SEPARATION JL 

(SNYDER, 1967, PAPER SHEETS, R<D.10) A 

(RUN II, 13 lb. SHEETS, R<0.10) O O 



A 
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TABLE III 



WEIGHT EFFECTS ON THE SLOPE 



SOURCE 








REJECTION FACTOR 






SLOPE 


Run 


I 


13 


lbs . 




0 


.05 








2.01 


Run 


I 


15 


lbs . 




0 


.05 








1.68 


Run 


I 


13 


lbs . 




0 


. 10 








1.62 


Run 


I 


15 


lbs . 




0 


. 10 








1.27 


Run 


II 


13 


lbs . 




0 


.05 








1.88 


Run 


II 


15 


lbs . 




0 


.05 








1.70 


Run 


II 


13 


lbs . 




0 


.10 








1.86 


Run 


II 


15 


lbs . 




0 


.10 








1.71 


. it 


seems to 


be 


a generally i 


accepted 


procedure 


to use 


the va 1 ue 


l. 10 


for the reject 


:ion factor. 


To s 


tudy 


the effec 


ts of 


c ha ng i ng 


re j( 


BCt 


ion 


factor , 


three values 


were 


considered : 


(i) R 


< 0.05, 


R < 0 


. 10 


, and 


(ii 


Li) R < 0.15. 


It 


was 


found that 


only 


the 



values of 0.05 and 0.10 would be useful. For both sheets A and 
B the values of k were increased by a factor of approximately 2.5 
when the rejection factor was increased from 0.05 to 0.10 (Table IV). 
Since the method of determining F(l) is only valid for R < < 1, this 
suggests that experimental values of k are high by at least a factor 
of two. 

The variation in the slope of the best-fit line, as a result of 
changing the rejection factor from 0.05 to 0.10, was not as con- 
sistent as that of k. The slope decreased in Run I but remained 
essentially constant in Run II (Table V). 
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TABLE IV 



EFFECTS 

SOURCE 


OF THE REJECTION FACTOR 

, t 2/3 -1 

k (cm sec ) 


ON k 

, , 2/3 -1 

k (cm sec ) 




(R < 0.05) 


(R < 0.10) 


Run I 


0.002 


0.004 


(13 lbs.) 


Run I 


0.001 


0.003 


(15 lbs.) 


Run II 


0.002 


0.005 


(13 lbs.) 


Run II 


0.002 


0.005 


(15 lbs.) 


EFFECTS OF 


TABLE V 

THE REJECTION FACTOR ON 


THE SLOPE 


SOURCE 


SLOPE 


SLOPE 




(R < 0.05) 


(R < 0.10) 


Run I 


2.01 


1.62 


(13 lbs.) 


Run I 


1.68 


1.27 


(15 lbs.) 


Run II 


1.88 


1.86 


(13 lbs.) 


Run II 


1.70 


1.71 


(15 lbs.) 



D # MEAN SEPARATION 

It was noted in the initial results that the separation between 
diffusers was not continuously increasing as was expected. Therefore, 
the mean separations of all pairs of diffusers were calculated for 
each pass of Runs I and II. The results of the mean separation 
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versus total time are shown in Figures 14 and 15. The variation of the 
mean separation with time suggested the possibility of wave action 
affecting the diffusion. To study the effects of- wave action the 
relationship between F(l) and wave spectra was determined. 



E. EFFECTS OF WAVE ACTION 



In this study let F (l) ~ where l is the separation between 

two particles, and the bar denotes an ensemble average. If the sea 
surface is disturbed by a monochromatic wave, then r)(x) = a sin (Hx-oot) 
where r\ is the free surface height, a_ is wave amplitude, H is wave 
number, x is horizontal distance, CD is angular frequency, and t is 
time. Consider two particles located on the x-axis and separated by 
a mean distance l. Then Iq = l + and l ^ = l + e^, where e is the 
change of separation due to wave motion. Since sampling time (minutes) 
is much larger than wave period (sec), and are independent. 

Thus, (At) 2 = (l 1 - l Q ) 2 = ( Sl - e 0 ) 2 



2 , 2 9 

“ e l + e 0 ‘ 2e l e 0 



= e l + e 0 * 



Thus, (Al) = = 2e since the mean quantities are independent of 

time. Here, the bar denotes either a time or an ensemble average. 

The value of e must now be calculated. Assume the particles to be 
located at mean positions 0 and l. As both particles move in a circu- 
lar orbit of radius <3, 

e(l, t) = a sin (-cot ) - a sin (h l - cot) 



and 

£ 2 2 2 

— = sin a) t + sin (nl-ot) + 2 sin cut sin (Kt-cDt). 
a 
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Averaging : 



€ | ,p - ^ - =1 + 2 sin cot sin (Hl-cot) . 
a 

Simplifying the equation using trigonometric identities gives 

- — ^ ^ = 1 - cos Hi. Consider l as some fraction (N) of the wavelength 
a 

(X): l = -NX . The equation now takes the form 



since X = Thus, 



£ ^ = 1 - cos 2 ttN 
a 



(i) 



2 2 

aT ■ 77 (! " cos 2 ttN). 



( 2 ) 



w At At 

Suppose that N is small. This is a reasonable assumption since a 

16-second swell has a wavelength of about 1200 feet. Then if l is 

approximately 200 feet, N~ 1/6, and 2ttN ~ 1. Since cos 1 = 0.54, 

2 

and the expansion cos x~ l~x /2 gives 0.5, we can use the expansion 
in (2), giving 

2 2 

9tt a 9 

(3) 






2 2 

which implies an N relationship (ie . Fal ). The power will decrease 
as N becomes larger (so that the approximation is no longer valid). 

The direction of the line between two particles must be considered. 
Normally this line will be at some angle (0) to the direction of wave 
travel. Then, 

e(l, t) = la sin (-cot) - a sin (Hi - cut)^ cos 8. 






When averaging over the direction, cos 0 = 1/2. Thus, (2) becomes 

2 

f w ~!af (1 - cos 2nN) > 



(4) 
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and (3) becomes 



‘ w At 



2 2 Q 

17 3 N 2 , N < 1/6. 



(5) 



Consider two wave trains at angles 0^ and 0^ to the line between 
particle^. Then, 

e(l, t) = a ^ /cn^, cos 0^ + a^ <jo) 2 , cos 0 2 . 



Since the waves are independent, 



e 2 = a^ 2 (1 - cos 2 ttN^) + a 2 2 (1 - cos 2rrN 2 ), 



( 6 ) 



Equation (6) can be generalized to any number of wave trains, M. Then 
using equation (5) gives 

2 M 9 9 

F (t) = fZa. N . 2 
w v ' At l l 



or 



2 M 

12 2 
F (l) = ~j~r~ Z h/ a. . 
w v ' 4At i l 



> 4 - 



(7) 



l ( 2 

In terms of energy spectrum S(h) 5 (7) becomes F^(l) = ~Kt ) K 
where the integral is the mean square surface slope. From (4), the 
exact form of (7) is 



M 

V l) = ^ S 3 i 2 (1 " COS K i l) ’ (8) 

Then the neighbor diffusivity due to wave motion can be calculated 
from a measured energy spectrum. 

The wave records at a position close to that where the experiments 

were conducted were available for the days the data were collected. A 

power spectrum was determined for each day using a digital computer 

(Figure 16). With the power spectra F (l) could be calculated using 

w 

various values of l (Appendix C). The power spectrum for Run II was 
chosen since it contained more energy than that for Run I. From 
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Figure 17, which represents the data for Run II in Appendix C, it was 
determined that the wave action had little direct effect on the 
diffusion of the particles. 
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FIGURE 16 
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APPENDIX A 

FLIGHT DATA 
Flight Data for Run I 
ALTITUDE (FEET) 

750 

875 

965 

1,100 

1,090 

1,180 

1,235 

920 

940 

1,045 

995 

895 

930 

1,000 

995 

1,000 

900 

955 

960 

975 



TIME INTERVAL (SEC) 

110 

112 

250 

112 

117 

126 

130 
511 
133 
136 
284 
165 
138 
122 
114 
127 
110 
259 

131 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 



Flight Data for Run II 
ALTITUDE (FEET) 

865 

815 

925 

865 

955 

1,060 

1,040 

900 

1,000 

940 

895 

860 

1,045 

920 

930 

950 

855 

860 

940 

940 

1,015 

1,010 

875 



TIME INTERVAL (SEC) 

110 

133 

122 

113 

128 

117 
122 

118 
119 

124 
132 
113 

125 
130 
105 
110 

98 

129 

128 

125 

135 

124 
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APPENDIX B 



RESULTS OF ANALYSIS 
Run I, 13 Pound Diffusers, R < 0.05 

2 

GLASS INTERVALS (m) NUMBER OF VALUES RESULTS (m /sec) 



65-80 


56 


F (71.3 1) 




0.01356 


80-90 


57 


F(85. 96)) 


= 


0.02444 


90-100 


44 


F (95 . 67 ) 


=r 


0.02500 


L00-LL0 


63 


F( 104 . 92 ) 


= 


0.03522 


110- 120 


52 


F(114.55) 


= 


0.03931 


120-130 


61 


F (125 .43 ) 


=r 


0.04186 


130-145 


64 


F(137 .53) 


= 


0.05098 


145-160 


45 


F( 15 1.82) 


= 


0.06979 


180-195 


46 


F ( 186 .49) 


= 


0. 10282 




Run I, 13 Pound Diffusers, 


R < 0..10 






40-55 


56 


F (45 . 87 ) 


= 


0.02329 


65-80 


85 


F (7 1. 87) 


= 


0.04522 


80-90 


75 


F(85 .81) 


= 


0.05129 


90-100 


53 


F (95.60) 


= 


0.05807 


100-110 


81 


F ( 105 . 07 ) 


= 


0.07598 


110-120 


69 


F (114 . 68) 


= 


0.09095 


120-130 


70 


F( 125 .35) 


= 


0.07648 


130-145 


90 


F(137 .51) 


= 


0. 12774 


145-160 


62 


F ( 152 . 11) 


= 


0 . 182 12 


160-175 


50 


F ( 165 . 58) 


= 


0.21077 


180-195 


56 


F ( 187 . 00) 


= 


0.18107 
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CLASS INTERVAL 


Run I , 


15 Pound 


Diffusers , 


R < 0.05 


(m 


2, \ 
/sec) 


(m) 


NUMBER 


OF VALUES 


RESULTS 


40-55 






59 


F (47 . 90) 


= 


0.00779 


55-65 






64 


F(60.52) 


= 


0.00884 


65-70 






56 


F(67 . 58) 


= 


0.01086 


70-80 






62 


F (74 .42) 


= 


0.01543 


80-90 






55 


F (85 . 84) 


= 


0.02056 


90-100 






53 


F (94.38) 


= 


0.02477 


100-115 






50 


F(107 .58) 


= 


0.02384 


115-130 






47 


F ( 122 . 59) 


= 


0.03485 


130-145 






56 


F( 136 . 97 ) 


= 


0.03984 




Run I , 


15 Pound 


Diffusers , 


R < 0.10 






25-40 






65 


F(34.22) 


=s 


0.01626 


40-55 






98 


F (47 . 90) 


= 


0.02275 


55-65 






95 


F(60.27 ) 


= 


0.02819 


65-70 






69 


F(67 .59) 


= 


0.02667 


70-80 






93 


F(74.31) 


= 


0.04279 


80-90 






79 


F (85 . 80) 


= 


0.05201 


90-100 






75 


F(94 . 47 ) 


= 


0.07250 


100-115 






69 


F ( 107 .43) 


= 


0.08183 


115-130 






55 


F(122 . 54) 


=S 


0.07943 


130-145 






64 


F(137 .01) 




0.06640 
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CLASS INTERVAL 
30-45 
45-60 
60-75 
75-85 

85- 95 
95-105 

105-115 

115-125 

125-135 

135-145 

145-155 

155-165 

165-175 

175-185 

185-200 

30-45 

45-60 

60-75 

75-85 

86- 95 
95-105 

105-115 

115-125 

125-135 



Run II, 13 Pound Diffusers, 


R < 0.05 






(m) NUMBER OF VALUES 


RESULTS 


(m ^/sec) 


39 


F(37 . 5 1) 


= 


0.00488 


42 


F(51.67) 


= 


0.01120 


45 


F(67 .50) 


= 


0.01312 


51 


F(80.33) 


= 


0.01982 


33 


F(90.09) 


= 


0.03118 


43 


F(99. 19) 


= 


0.02912 


59 


F(109.85) 


= 


0.04057 


53 


F(120 . 11) 


= 


0.02847 


61 


F(130.36) 


= 


0.03986 


55 


F(139.87) 


= 


0.05761 


67 


F(150 .34) 


= 


0.07082 


58 


F(160 . 84) 


= 


0.08120 


40 


F ( 169 . 51) 


= 


0. 12072 


43 


F(180.41) 


= 


0.09218 


41 


F(191.28) 


= 


0.10879 


Run II, 13 Pound Diffusers, 


R < 0.10 






81 


F(37 .62) 


= 


0.01185 


72 


F(52 . 10) 


= 


0.03320 


83 


F(67 .81) 


= 


0.05298 


78 


F(80.33) 


= 


0.05729 


55 


F(89 . 86) 


= 


0.09132 


70 


F (99 . 69) 


= 


0.10325 


91 


F ( 109 . 64) 


= 


0.11272 


74 


F(120.06) 


= 


0.10330 


95 


F(130.23) 


= 


0.15158 
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CLASS INTERVAL 


(m) 


NUMBER OF VALUES 


RESULTS 


(m ^/sec) 


135-145 




88 


F(140.08) 


= 


0.20152 


145-155 




97 


F(150.09) 


= 


0.21462 


155-165 




91 


F(160.26) 


= 


0.24438 


165-175 




64 


F(169.53) 


= 


0.32294 


175-185 




54 


F ( 180.44) 


= 


0.20824 


185-200 




58 


F ( 191.42) 


= 


0.28200 




Run II, 


15 Pound Diffusers, 


R < 0.05 






40-55 




45 


F(48. 54) 


= 


0.00967 


55-65 




52 


F (59 . 90) 


= 


0.01271 


65-75 




59 


F(70.01) 


= 


0.01323 


75-90 




53 


F(80.81) 


= 


0.02023 


90-105 




60 


F(98. 94) 


= 


0.03340 


105-110 




39 


F ( 107 . 53) 


= 


0.04644 


110-115 




48 


F(112.29) 


= 


0.04404 


115-120 




45 


F (117 .78) 


= 


0.03473 


120-125 




38 


F ( 122 .73) 


= 


0.04667 


125-130 




50 


F(127 . 10) 


= 


0.05341 


130-140 




59 


F(134.89) 


= 


0.04745 


140-150 




55 


F ( 144. 84) 


= 


0.05272 


150-160 




62 


F(155.33) 


= 


0.07673 


160-170 




39 


F(164. 89) 


= 


0.06392 


170-185 




57 


F ( 177 .51) 


= 


0.09138 


185-200 




38 


F(192 .15) 


= 


0.07876 
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CLASS INTERVAL 


Run II , 


15 Pound Diffusers, 


R < 0.10 


(m ^/sec) 


(m) 


NUMBER OF VALUES 


RESULTS 


40-55 




87 


F(48.43) 


= 


0.03174 


55-65 




87 


F (59 . 65) 


= 


0.04248 


65-75 




92 


F(70. 14) 


= 


0.05138 


75-90 




80 


F(80.60) 


= 


0.06454 


90-105 




89 


F (98 . 57 ) 


= 


0.08951 


105-110 




57 


F(107 .58) 


= 


0.10190 


110-115 




79 


F(112 .48) 


= 


0.14685 


115-120 




67 


F(117 .65) 


= 


0.12690 


120-125 




64 


F(122 .67) 


= 


0.14780 


125-130 




70 


F(127 . 19) 


= 


0.15046 


130-140 




85 


F(135.11) 


= 


0.16463 


140-150 




79 


F(145. 11) 


= 


0.17404 


150-160 




86 


F(155.20) 


= 


0.18276 


160-170 




64 


F(164.92) 


= 


0.29894 


170-185 




83 


F(176.35) 


= 


0.26843 


185-200 




51 


F ( 192 . 13) 


= 


0.30326 
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APPENDIX C 



t (m) 


RUN II WAVE SPECTRUM RESULTS 

F ( t ) (m 2 /sec) 
(x 10' 4 ) 


10.0 


0.00617 


15.0 


0.01290 


20.0 


0.02085 


30.0 


0.03741 


40.0 


0.05350 


50.0 


0.07032 


60.0 


0.08748 


70.0 


0.10236 


80.0 


0.11364 


90.0 


0. 12175 


100.0 


0.12645 


110.0 


0.12682 


120.0 


0.12356 


150.0 


0.10919 


170.0 


0.09773 


190.0 


0.08814 


210.0 


0.08369 
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